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Abstract 

We study the renormalization of the AS" = 1 effective weak Hamiltonian with overlap 
fermions. The mixing coefRcients among dimension-six operators are computed at one loop 
in perturbation theory. As a consequence of the chiral symmetry at finite lattice spacing and 
of the GIM mechanism, which is quadratic in the masses, the K tttt and if — > tt matrix 
elements relevant for the A/ =1/2 rule can be computed without any power subtractions. The 
analogous amplitudes for e'/e require one divergent subtraction only, which can be performed 
non-perturbatively using Jf — > matrix elements. 
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1 Introduction 



The non-leptonic AS = 1 weak transitions are far from being understood. It is well known that the 
A/ = 1/2 amplitudes^ are enhanced (the so-called A/ = 1/2 rule) with respect to the corresponding 
AI = 3/2 ones (by roughly a factor 20 in ^ vrvr decays); the latest measurements of e'/e [Q] 
confirm the rather large value found by NA31 Q, and the up-to-date world average is roughly 
2 • 10~^. Given the present limited knowledge of the long-distance QCD effects in non-leptonic 
AS* = 1 amplitudes, these experimental results cannot be predicted from first principles and well- 
controlled approximations. 

The short-distance QCD corrections to the AS* = 1 effective Hamiltonian can be reliably 
computed in perturbation theory^ p. In the Standard Model (SM) the Wilson coefficients are 
known up to the Next-to-Leading Order (NLO) 0]. They account for a factor two of the 
enhancement in the AI = 1/2 rule and cannot reproduce the experimental value of e'/e, when 
ImAt = ytdYts is determined from the unitarity triangle and the Vacuum Saturation Approximation 
(VSA) estimate of {ttttIOgIK) (see below) is used |]-|0|. The SM could expl ain both effects if 



penguin contractions, neglected in the VSA, give contributions to the relevant matrix elements 
definitely larger than their factorized values . Up to now no reliable first-principle computations 
of the AI = 1/2 matrix elements exist. 

Lattice QCD is the only known method which allows one to compute non-perturbative QCD 
contributions to physical amplitudes from first principles. In the most popular lattice regulariza- 
tions, i.e. Wilson and staggered fermions, the main theoretical aspects of the renormalization of 
composite operators are fully understood [11|-|13]. The results of quenched calculations of the ma- 
trix elements relevant for the K^-K^ and B^-B^ mixings have reached a high level of accuracy and 
are commonly used in phenomenological analyses of the unitarity triangle |14]-[^]. These lattice 
results have been crucial in constraining the parameters of the CKM matrix in the Standard Model 
(SM) [l^, P and beyond |l^. Techniques to compute AI = 1/2 amplitudes for the AS = 1 Hamil- 



tonian have been developed for both Wilson and staggered fermions |21[], but these methods 
have not yet produced useful results^. There are two major difficulties: 

• In Euclidean space there is apparently no simple relation between the two- (many-)body tran- 
sition matrix elements and the corresponding correlation functions at large time separation 
[23|. This problem is common to all lattice discretizations and there are proposals to solve it 



either by using finite volume techniques or studying the correlation functions in the large 
volume limit p6| ; 

• The operators we are interested in mix with lower-dimensional operators with coefficients 
which can diverge as inverse powers of the lattice spacing. These mixings can be more severe 
when chiral symmetry is broken by the regularization, like for Wilson fermions. 



^ In this paper we will study only amplitudes for K tttt decays. 



■^In the last few years there has been much progress in constructing a non-perturbative formulation of chiral gauge 
theories on the lattice ^, ^. 

■^Attempts with domain- wall fermions are in progress ]22|]. 
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In the last few years it has been understood [27|-|29] that chiral and flavor symmetries can be pre- 
served simultaneously on the lattice, without fermion doubling, if the fermionic operator D satisfies 
the Ginsparg- Wilson Relation (GWR) []30|| . This implies an exact invariance of the fermionic action 
which can be interpreted as a lattice realization of the standard chiral symmetry at finite cutoff 
[p^j . However the GWR itself does not guarantee locality, the absence of doubler modes and the 
correct classical continuum limit. 



Neuberger, through the overlap formalism |^], found a solution [27| of the GWR which satisfies 
all the above requirements and is localQ |^2|. The complicated form of Neuberger 's operator renders 
its numerical implementation quite demanding for the present generation of computers. However, 



some progress has been achieved [33] and Monte Carlo simulations are already feasible, at least in 
the quenched approximation. 

In this paper we study the renormalization of the AS = 1 effective weak Hamiltonian in the 
overlap regularization. Neuberger's action implies many theoretical advantages: 



• The quark masses renormalize only multiplicatively, the mixings among operators with dif- 
ferent chirality are forbidden p^, |33] and their parity-violating and parity-conserving com- 



ponents renormalize in the same way |34]; 



• The GIM mechanism is as powerful as in the continuum to suppress mixing with lower- 
dimensional operators; 

• The action and therefore the spectrum of the theory are free from 0(a) discretization errors. 



The 0(a) improvement of the local fermionic operators is greatly simplified |35|. 



As a consequence the K — > vrvr and, if chiral perturbation theory (xPT) is used, X — > vr matrix 
elements for the A/ = 1/2 rule can be computed without any power subtractions. The anal- 
ogous matrix elements for e'/e require only one divergent subtraction which can be performed 
non-perturbatively using on-shell K ^ matrix elements. On the contrary, when the regulariza- 
tion breaks chiral symmetry the K ^ ir matrix elements for the AI = 1/2 rule and for e'/e require 
one and two power-divergent subtractions respectively. 

We also compute the mixing coefficients among dimension-six operators at one loop in per- 
turbation theory^ A further determination of these renormalization constants can be obtained 
using numerical non-perturbative methods such as those in Refs. |3^. However a perturbative 
analysis is useful for analytically studying the mixing pattern of the operators, often gives accurate 
estimates of the renormalization constants and furnishes a consistency check of the non-perturbative 
methods. 

The paper is organized as follows: in section |^ we define the A5 = 1 effective Hamiltonian in 
the continuum; in sections |3| and ^ the overlap fermion action and the lattice A5 = 1 operators are 

The overlap-Dirac operator is not ultra-local. The Neuberger kernel satisfies a more general definition of locality, 
i.e. it is exponentially suppressed at large distances with a decay rate proportional to 1/a. 

^These renormalization constants for domain-wall fermions in the limit of infinite fifth dimension have been 
computed in j36j. 
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introduced, and in section |5| we compute at one loop the penguin diagrams of these operators. In 
sections ^ and ^ we describe the renormahzation pattern of the relevant operators for the A/ = 1/2 
rule and e'/e and we report the main results of this paper; in section ^ we state our conclusions. 



2 The AS = 1 effective Hamiltonian in the continuum 

The AS = 1 effective Hamiltonian above the charm threshold is given by 

n^ff' = A„^[(l-r)(Ci(/i)(Qi(/x)-Q^(/x))+C2(/i)(Q2(/i)-Q§(/^) 

10 

+ Tj2a{fs)Q,{fi)] , (1) 

i=l 

where Gp is the Fermi coupling constant, r = —Xt/Xu, Xu = VudV*g, Xt = VtdV^* and fj, is the 
renormahzation scale of the composite operators. The Wilson coefficients Ci{^) at the NLO in 
QCD and QED have been given in |^ . The operator basis is 

Q2 = {sj1u){uj1d) 

Q 

3 

^ q 

^ q 

Q\ = (s-"7^c^)(c^^^i") 

= (.-7^)(c-7^d) , (2) 

where ^ = 7'^(1 =F Ts); are the quark electric charges and a, (5 denote color indices (which are 
omitted for the color-singlet operators). Below the bottom threshold, q runs over the active flavors 
u, d, s, c. Other operators would appear in Eq. (||), but in the Standard Model they are numerically 
irrelevant and they will not be considered in this papei]^. 

The electro-penguin operators Qi {i = 7, . . . , 10) give significant contributions to e'/e only 
through their A/ = 3/2 components. In the limit = m^, do not mix with lower-dimensional 
operator^ and their renormahzation constants in the overlap regularization can be obtained from 

®In some supersymmetric extensions of the SM, large contributions to e'/e can be induced by the magnetic and 
chromomagnetic operators [ ^ . 

^Estimates for the matrix elements of these operators in the quenched approximation have been obtained on the 
lattice in the standard regularizations |14[- jl7[. 



[^]. Therefore in the fohowing we will focus only on the computations of the matrix elements of 
Qi, . . . ,Q6 (and Qi.Qi)- 

Under flavor SU{3)l x SU{3)f(, transformations, the combination 

Q27 = 2Q2 + 3Qi -Q3 + QI (3) 

transforms as a {27l,1f{) multiplet and induces both AI = 1/2 and AI = 3/2 transitions. The 
operators 

Q- = Qi- Q2 (4) 

and Qi {i = 3,4,5,6) transform as {SlAr) operators and induce pure AI = 1/2 transitions. 

The octet operators are linearly dependent; by performing Fierz transformations it is easy to 
show that 

Qa = Qz + Q2-Qi + Q'2-QI. (5) 

For ^ < rric the operator basis is given by Qi — Qio and the sums over q in the penguin operators 
run over n, d, s only. The relation (|5[) becomes 

Q4 = Q3 + Q2 - Qi (6) 

and in Eq. Q\ has to be dropped. 



3 Basic definitions for overlap fermions 



The QCD lattice action we consider for massive fermions is 
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1 - -TV 
6 



Up + uj. 



+ 



1 - —amgjDiy + rUg 



(7) 



where, in standard notation. Up is the Wilson plaquette, go is the bare coupling constant and rUg 
is the bare fermion mass. is the Neuberger-Dirac operator defined as 

1 



X 



a V VXTX 
Dw P , 



a 

1 r 



(8) 



where 

Dw = §7,.(V^ + V*) - -aV;V^ (9) 

is the Wilson-Dirac operator^, < r < 1 and < p < 2r. and V* are the forward and backward 
lattice derivatives, i.e. 

U^{x)q{x + afi) - q{x) 



a 
1 r 



- UUx - afi)q{x - afi) 



(10) 



All perturbative computations reported in this paper are performed with r — 1. 
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where U^{x) are the lattice gauge hnks. 

The overlap-Dirac operator satisfies the Ginsparg- Wilson relation 

a 

'J^Dn + DnJ5 = -Dn^^Dn (11) 

p 

which, in the massless limit, implies a continuous symmetry of the action which may be interpreted 
as a lattice form of the chiral invariance |29|: 

= %q , = qj5 , (12) 

where 75 is defined as 



and satisfies 



If we define the chiral projectors as 



75 = 75(1 --^jv) (13) 
7I = 75 , 7i = 1 • (14) 



P± = ^(l±75), P± = ^(l±75), (15) 
the fermionic fields can be decomposed into left- and right-handed components 41 1 

gR,L = P±q qR,L = qP^ (16) 

which transform under lattice chiral rotations in the same way as the corresponding fields in the 
continuum. The massive fermion matrix in Eq. ^ can be written as 

DN{m) = P+DnP- + P-DnP+ + m{P_P_ + P+P+) . (17) 

The composite fermionic operators can be defined to have the same combination of left and right 
components as in the continuum. This guarantees that under lattice chiral rotations they transform 
as the corresponding operators in the continuum. As an example, the bilinear operators can be 
defined as 

Or = gif (?2 , (18) 



where T = T{1 — ^Dn)- The on-shell matrix elements of these operators are 0{a) improved ||3J 
The Feynman rules corresponding to the action in Eq. are given in |34] and are in agreement 
with those used in |p^ -[44|. Reference |34] also gives in detail all conventions and symbols used in 
this paper. 



4 = 1 four-fermion operators on the lattice 



The four-flavor QCD lattice action in the overlap regularization reads 



S 



-y 



1 - -Tr 
6 



Up + ul 



+ 



i:4( 



1 - -^arUq^DN + rriq 



(19) 



where q = u,d, s, c. In the hmit = m^i = rUg = 0, this action is invariant under the SU{3)l x 
SU{3)ii group and the corresponding transformations are defined as 

ipR ^ V4^R , i^R ^ i^RV^ , (20) 

where ipL,R = {ul^r, di^R, sl,r) is a cokimn vector in flavor space and 

C/t[/ = y^v = 1 

detf/ = dety = 1 . (21) 
When u,d,s acquire a mass, the action is invariant if the hght quark masses are transformed as 



/ m 



M 



\ 



u 

rud 



^ [/mW . (22) 
A convenient definition for the lattice four-fermion operators of the AS" = 1 effective Hamiltonian 



IS 



Qi = (s-"7^^)(n'^7£c^") 

23,5 = {si1d)Y,{qjlj,q) (23) 
q 

Q4,6 = {-s"rLd^)T.('l^llRqn 

g 

Ql = (.-7^c^)(c^7^d") 
= {slicmid), 



where 7^ ^ = Ilr^^ ~ 'Tp^^)- Under the rotations in Eq. (20), these operators transform as the 
analogous operators in the continuum and their on-shell matrix elements are 0{a) improved. 

Under renormalization the operators in Eqs. ( ]23| ) can mix with equal- or lower-dimensional 
operators which have the same transformation properties under SU{3)l x SU{3)r or, due to the 
explicit soft chiral symmetry-breaking term, with multiplets of different chirality properly multiplied 
by light quark masses. As in the continuum, Q27 is multiplicatively renormalizable while the 
(8l, 1_r) operators can mix among themselves as well as with two lower-dimensional operators 
[ p] , p5[ | (besides operators which vanish by the equations of motion) 

Qcr = 50 {ms + md)sa^yF^yd+{'ms-md)sa^uF^^d (24) 
Qm = (ms + md)sld + {nis - md)s%d , (25) 

where F^i, and F^i, are the non-abelian gluon field tensor and its dual, a^^ = — ^Dn), 

75 = 75(1 — ^^-^n) and 1 = (1 — ^F)^). The definition of the properly subtracted operators 
necessary for the AI = 1/2 rule and for e'/e will be given in sections ^ and |^. 
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5 Penguin diagrams in perturbation theory 



In this section we discuss the calculation of the penguin diagrams at one loop in perturbation theory. 
The four-point Green's functions of a four-fermion operator between quark states are defined as 

Gi{xi,X2,Xz,Xi) = {qi{xi)q2{x2)Qi{^)q?,{x:i)qi{x4)) , (26) 

where xi, 3:3 and X2, X4, are the coordinates of the outgoing and incoming quarks respectively. Note 
that Gi depends implicitly on the color and Dirac indices carried by the external fermion fields. 
The Fourier transform of the Green's functions (^), with external momenta p and p' chosen as in 
Figs. |l] and is defined as 

Gi{p,p') = j dxidx2dx3dx4e-'P^'^'-^*^e-'P'^^^-'^^'^Giixi,X2,X3,X4) ■ (27) 

The corresponding amputated correlation functions are defined as 

Up,p') = S'\p)S'Hp')G,ip,p')S-\p')S-Hp) , (28) 

where S{p) denotes the quark propagator in momentum space. At one loop, the only diagrams 
which contribute to the mixing among dimension-six operators are the penguin diagrams in Fig. |^ 



and the vertex diagrams in Fig. As for the AF = 2 operators in Ref. |34|, it can be proven that 
the operators in Eqs. ( |23| ) renormalize as the corresponding local ones. Therefore in the following 
all one-loop computations are performed with local operators. 

The vertex diagrams in Fig. |^ have been computed in Ref. |3^; they are sufficient for the 
renormalization of the AF = 2 operators and do not give rise to mixing with lower-dimensional 
operators. The penguin diagrams in Fig. |l| are computed here for the first time and constitute one 
of the main results of this paper. The calculations involve the gauge-invariant quantity 



Figure 1: Penguin diagrams for four-fermion operators. 




V,{ak,am) = (^4°^^ " + |' ^ " + ^) > k = p' - p , (29) 
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Table 1: Perturbative values of the finite parts of the amputated correlation functions in the Landau 
gauge for the self energy, the scalar density and the vector current. The gauge-invariant penguin 
constant -Bpeng defined in Eq. ( |30[ ) is reported in the last column. 



which does not depend on the flavor, spin and color structures of the operators, and is the result 
of the internal quark loop once these structures have been factored out. The tree-level overlap 
propagator 5^ 



•^'^\p) and the quark-gluon vertex vli^\p,p') 



' ! ^'^^ deflned in Ref. ||3^. Eventually we 

will impose the RI renormalization scheme and therefore we can perform all the computations 
in the massless limit. V^{ak,0) has to be computed by expanding vertices and propagators to 
second order in ak, since in the diagrams in Fig. || it is contracted with the gluon propagator 
G{ak) ~ (ak)-^. 

For the analytic computations we have used FORM codes which we have developed specifically 
for this problem; the output has been integrated on a 60^ grid. Among other checks, we have 
computed Vfj,{ak,0) using different assignments for the loop momenta. We also checked that the 
spectator quark line contributes only with its continuum expression. The result for overlap fermions 
is 



V^,{ak,0)=i 



go 

167r2 



1 



loga^/c^ + B 



peng 



(30) 



where Bpeng is given in Table |l] for several values of p. Notice the structure in k which comes from 
gauge invariance. The ^ term vanishes, up to 0(5q), if the equations of motion are used for the 
spectator quarks |4(:], and the k'^ in the remaining structure 'y^k'^ disappears when combined with 



the gluon propagator. 
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We have used our codes to compute the penguin diagrams also in the Wilson case, 



: go 

'I67r2 

: go 

'I67r2 



J log a^A;^ + 0.5709352 



1 , a^P 



0.1922181 



(31) 



and we are in agreement with the results by Bernard et al. [^]. It is interesting to note that 
the overlap result for = 1 is not too far from the Wilson number. In the Wilson case the lack 
of chiral symmetry causes, even in the massless case, the presence in V^{ak,0) of an additional 
power-divergent term which contains tensor operators: 

- 0.39847257 i^a^A^A- (32) 

We have also observed that, as expected, for overlap fermions the coefficient of the tensor term is 
zero (within our integration precision). This constitutes another test of the numerical accuracy of 
our integrals. 

Figure 2: Feynman graphs for the vertex corrections to four-fermion operators. 





A generic local four-fermion operator can give rise to one or more of the following penguin 
types^: 

Pai = {srrq){qT'r'd) 

^We use the notation sVTd = sfVijTaijdP- , V and V are general Dirac matrices, T and T' axe color matrices (1 
or f^), and q, (3 are color indices. 
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Pa2 = (^«rg/^)(g^r'd") 

Pbi = {sTTd){qT'T'q) (33) 

Pb2 = (s°rd^)(g^r'(?°). 

The Wick contractions which form the quark loop are always meant to be only between q and q, 
and at one loop they give 



lAi = -igoG{ak){sT V^{ak,am)T' rt^T d){ft^-f^f) 

Ia2 = -i9oG{ak){sr V^,iak,am)r' tr(t^) d){ft^j^f) 

TBi = -igoG{ak){sT ti{-V^{ak,am)T') Ttr(t^T') d){ft^j^f) 

Ib2 = -igoG{ak){sT tv{-V^{ak,am)T') d)ift^-/^f) 



(34) 



where G{ak) is the gluon propagator given in Ref. and / is the spectator quark. For the 
contractions Pai and Pa2 tbe relevant diagram is the one on the right of Fig. |l], while for Pbi 
and Pb2 it is the diagram on the left. Weak operators in which q = d oi q = s, since they have 
three quarks of the same flavor, can undergo two Wick contractions instead of only one, therefore 
contributing with both diagrams to their renormalization factors. 

The 1-loop penguin corresponding to Za2 is zero because of the tracelessness of the matrices; 
when T = T' = 1 the penguin contribution 2bi also vanishes for the same reason. Notice also 
that, since Vf^ enters in the expressions (^) either as TVpT' or as rtr(— "P^F'), the contribution of 
the CT^A term of Eq. (|3^ ) in the Wilson case vanishes for all operators considered in this paper. 



6 Renormalization pattern for the AI = 1/2 rule 

For renormalization scales above the charm mass the bulk of the CP-conserving AI =1/2 ampli- 
tudes comes from the operators 

0± = (Qi - Ql) ± (Q2 - Q'2) ■ (35) 

The contributions which arise when the top quark is integrated out are suppressed by a factor 
Xt/Xu and can be neglected[^ 

Chiral symmetry forbids mixings with other dimension-six operators and the flavor structure 
forbids mixings between Oj[- and 0-. When chiral symmetry is preserved, the GIM mechanism 
ensures that the mixing with Q^- and Qm (Eqs. ( p4[ ) and (^)) vanishes as —rn^ when rric —>■ ruu- 
This can be shown as in the continuum: the penguin contractions of left-left four-fermion operators 
take non-vanishing contributions only when the helicity, which is preserved through the massless 
quark lines, is flipped by an even number of mass insertions. Therefore in the overlap regularization 
the renormalized operators read as 

0±(Ai) = Z±{^^a,g^o)^±{a) + O{a^) (36) 

0±{a) = ±{a) + {ml - ml)C^{gl)QM , 

^°The case with a charm integrated out can be easily obtained from the discussion in the next section. 
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where Z±{fj,a, g^) are 0(1) logarithmic divergent renormahzation constants, C™(5q) are finite mix- 
ing coefficients which are suppressed by a factor Ug and the contributions proportional to 



Qo- give 0{a ) effects. Z± and C™ can be determined by imposing two renormahzation condi- 
tions which, in perturbation theory, would specify their two- and four-quark amputated Green's 
functions. 

For Z± we opt for the RI renormahzation scheme proposed in |37, 48 1. Z^^ are defined by 
imposing that the renormalized amputated Green's functions at large fixed Euclidean scale in 
the Landau gauge and in the massless limit, are equal to their tree- level values^, i.e. 



(^a,go')A±(p,p'; 



(37) 



At one loop the relevant contributions to Z^^ come from the vertex diagram in Fig. H while the 
penguin diagrams cancel in the massless limit due to the GIM mechanism. Using the results in 
Ref. \M] we obtain 



Z^' 



where 



7- 



167r2 



(0) 

(0) _ 



(0) 

7± 



log(/x^a^) + B± + 2 



Ni 



1 



-B,i 



(38) 



B. 
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By - i:Bs 



-XBs + Bv) . 



(39) 



The constants B^, Bs and By, which enter the renormahzation of the quark propagator, the scalar 



density and the vector current, were computed in |43, 34| for several values of p and are reported in 
Table |l| for completeness. The anomalous dimensions 7^^*^ are in agreement with those in Refs. ||6|, ^ 
and the matching coefficients between RI and one of the MS schemes can be found in |48|. 



The mixing coefficients Cm^ can be determined by imposing suitable renormahzation condi- 
tions on two-quark correlation functions computed at two loops in perturbation theory. As in 
the continuum, these coefficients are not needed for the physical K — > vrvr matrix elements. For 
— > TT amplitudes we will show that they can be fixed non-perturbatively by imposing suitable 
renormahzation conditions on the K ^ matrix elements. 



The K — > vrvr matrix elements 

{n7T\d±{p)\K) = Z±(/.a,5o')(vr^| fo±(a) + (ml - ml)C^{gl)Q. 



\K) 



(40) 



can be extracted directly from the four-point functions once the final state interactions have been 
properly taken into account |24]-|26|. Following the on-shell matrix elements of the pseu- 
doscalar density V = s%d satisfies the Axial Ward Identity 

ZA{TnT\d^Af,\K ) = (m, + md){TnT\V\K) + 0{a^) , (41) 

^'^In the Landau gauge this is equivalent to the use of the projectors as in Ref. | ^^ . 

^^Throughout this paper the Ai(p,p')'s for the four-fermion operators include only diagrams which contribute to 
the mixing among dimension-six operators. 



11 



where Af_i = s'yfj,%d and Za is its renormalization constant, and up to O(a^) effects {ttttIVIK) 
vanishes when the momentum Ap inserted by the operator is zero. Eventuahy the physical matrix 
elements in the RI scheme read 

{7T7r\d±if,)\K) = Zfifia,gl){7T7r\0±ia)\K) + ©(g^) . (42) 

Even if the momentum inserted is not zero, the contribution proportional to Ap is finite and 
suppressed by a factor a^. Therefore the extrapolation to the physical point Ap = would not be 
a problem. The direct computation described above does not rely on chiral perturbation theory 
and in principle the use of 0{a) improved operators renders the computation of the matrix elements 
sensitive to higher-order chiral contributions. The two main advantages of the overlap regularization 
with respect to the Wilson fermions are that the mixings with dimension-three operators are finite 
even at non-zero momentum transfer (Ap ^ 0) and the improvement of the bilinear and four-fermion 
operators can be performed without tuning any parameters. The main difficulties of this method 
are that the final state interactions have to be taken properly into account and the computations 
of the four-point functions are technically quite complicated which could lead to noisy signals. 



These difficulties can be avoided if one applies the method proposed in pO| , 21 1, where it has 
been suggested to use chiral perturbation theory to relate K vrvr to ^ vr matrix elements and 
to compute the latter on the lattice. At the leading order in ^PTf^ 



(7r+(p)|0±|i^+(g)) = a±p-q-6±^, (44) 



{0\O±\K^) = ib^(m\-inl) (43) 

U 

(7r+7r-|O±|if0) = ia^ ""^"""' , (45) 

Jtt 

and a± can be computed by extracting {tt^ {p)\0±\K^ (q)) from the three-point functions and 
fitting the results as a function oi p ■ q and m?^. The unphysical contribution proportional to 5± 
is not divergent. This is one of the main results of this paper. 5± can also be subtracted by fixing 



in such a way that |45| 



(0|O±|k0)=0, (46) 



I.e. 

Op{mj. - mi) rus - rud 

with 

{0\V\K^) = i6p . (48) 

The comparison of the two methods is a check of the whole procedure. Since only single-particle 
states are involved there are no problems with final state interactions and the three-point functions 
are technically easier to simulate. 

From the considerations reported above, we can appreciate the advantages of using Neuberger's 
fermions to compute the matrix elements necessary for the AI = 1/2 rule: 

^^For higher order results see Ref.fcol]. 
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• the GIM mechanism combined with chiral symmetry constrains the mixing coefficients with 
lower-dimensional operators to be proportional to — m^; 

• the mixing coefficients with parity-conserving and parity-violating components of lower- 
dimensional operators are the same and are multiplied by — and + md respectively; 

• 0=^(0) do not mix with dimension-six operators with different chiralities; 

• the 0{a) improvement of the operators does not require any tuning of parameters. 

On the contrary, in the Wilson regularization the parity-conserving components of can mix with 
dimension-six operators of different chiralities, their mixing with lower-dimensional operators is no 
longer proportional to -|- and the GIM mechanism gives only a factor nic — rriu- Therefore 
the contributions of the chromomagnetic operator have to be subtracted and the scalar operator 
mixes with quadratically divergent coefficients. Finally for Wilson fermions the 0(a) improvement 
of the four-fermion composite operators requires tuning many parameters. 



7 Renormalization pattern for e'/e 

For CP violating processes in kaon decays, when the top quark is integrated out, the GIM mech- 
anism is not operative and the presence or not of an active charm does not substantially modify 
the mixing pattern of the relevant operators. Therefore, in this section we will consider the case 
with the charm quark integrated outp^. The corresponding AS" = 1 effective Hamiltonian can be 
obtained from Eq. (|l]) and is reported, for example, in Ref. |^. 

Analogously to the continuum, the operator 

Q27 = 2Q2 + 3Qi - Q3 (49) 
does not mix with lower-dimensional operators and it is multiplicatively renormalizable: 

Q27(/^) = 2'27(/ua,c/o) 227(0) • (50) 

At one loop in perturbation theory only the diagrams in Fig. |2| give contributions to Z27{na, g^) 
and in the RI scheme, defined analogously to Eq. (p7|), we obtain 



167r2 



^log(/.V)+i?27 + 2^5^ 



(51) 



with 



(0) 
727 



B27 = -^By - -Bs 



(52) 



The case with a dynamical charm can be easily obtained from the results reported in this and in the previous 
section. 
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Also in this case chiral perturbation theory can be used to relate K vrvr to K ^ it matrix 
elements and the latter can be extracted from the three-point functions computed numerically. 

The mixing pattern of the QCD penguin operators 



23.5 = (s7L'^)Z!W,i?9) 

24.6 = {s^l1d(')Y,{q^-fl^q^ 



(53) 



IS 



Qj{a) 



Zij{na,gQ)Qj{a) 



z,j = 3, ... ,6 



(54) 



Q,(a) + CJigDQAa) + C™(5o')Q™(a) . 



m/2\ 



At one loop the overall renormalization constants Zij{^a, Qq) take contributions from the penguin 
and vertex diagrams in Figs. || and |2[ Imposing the RI renormalization conditions 



we obtain 



Zf^{lxa,gl)K,{p,p') 



167r2 



with i, j = 3, . . . 6, 



;{o) 



/ 22 

. 2/ 



2/ 



^\og{^?a^)+B,,+25i, ' 



22 

y 

-2 + 



2/ 





2/ 
3 



4 
~9 
2/ 
' 9 

2 

2/ 



2iV, 



4 
3 

2/ 
3 



-6 



-16 



2/ 



(55) 
(56) 



(57) 



ll 
3 



5 



3 



By - Bs + —Bpeng 







V 



3 -Openg 



-By — Bs — 2B. 



peng 



g-By + - f Bpeng 



peng 



2 ^peng 

Lb 

3 -Cpeng 



-2B 



peng 



f Bpeng 
—By + -B5 

- f Bpeng j 



(58) 



where[^ / = 3. Notice that all the mixing coefficients in Eq. ( [5^ ) depend on four constants only. 
^^The case of a dynamical charm is obtained with / = 4. 
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The mixing coefficients C- are finite, can be computed in perturbation theory and are relevant 
only if one is interested in higher-order contributions in xPT. We postpone the computation of C" 
to a future publication. 

The mixing coefficients CJ" are quadratically divergent in the lattice spacing a and have to 
be computed non-perturbatively. As in the previous section, they can be fixed by imposing the 
renormalization condition p5[ 

(OlQ.li^O) =0. (59) 

The K — > vrvr or, if xPT is used, the i^T — > vr matrix elements can be computed analogously to 
section ^. The K ^ n matrix elements require one power-divergent subtraction only and there 
are no extra mixings with dimension-six operators of different chirality. For Wilson fermions, on 
the contrary, mixings with dimension-six operators of different chiral behavior are allowed and two 
power-divergent subtractions are required. 

8 Conclusions 

We have analyzed various possibilities for computing the matrix elements relevant for the AI = 1/2 
rule and for e'/e with Neuberger's fermions. When chiral symmetry is preserved at finite lattice 
spacing the renormalization factors of the parity-violating and parity-conserving components of 
the four-fermion operators are the same. Therefore, contrary to the Wilson case, there are no 
extra difficulties in extracting the matrix elements of the properly renormalized parity-conserving 
operators. We have shown that the K tttt and, if chiral perturbation theory is used, X — > tt 
matrix elements for the AI = 1/2 rule can be computed without any power subtractions. This is 
a consequence of the exact chiral symmetry at finite lattice spacing and of the GIM mechanism 
which is quadratic in the quark masses. The analogous matrix elements for e'/e require only a power 
subtraction which can be performed non-perturbatively using on-shell K ^ matrix elements. All 
the required renormalization coefficients among dimension-six operators have been computed at 
one loop in perturbation theory and depend on four numerical constants only. We believe that the 
overlap-Dirac operator is a promising regularization to solve these long-standing problems. 
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